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Abstract. A hexagonal patch is a plane graph in which inner faces have length 6, 
inner vertices have degree 3, and boundary vertices have degree 2 or 3. We con- 
sider the following counting problem: given a sequence of twos and threes, how 
many hexagonal patches exist with this degree sequence along the outer face? 
This problem is motivated by the enumeration of benzenoid hydrocarbons and 
fullerenes in computational chemistry. We give the first polynomial time algo- 
rithm for this problem. We show that it can be reduced to counting maximum 
independent sets in circle graphs, and give a simple and fast algorithm for this 
problem. 

Keywords: graph algorithms, computational complexity, counting problem, pla- 
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1 Introduction 

The notions used and problems introduced in this section are defined more formally 
in Section 2. A plane graph G is a graph together with a fixed planar embedding in 
the plane. The unbounded face is called the outer face and the other faces are called 
inner faces. The boundary of the outer face is simply called the boundary of G. A 
hexagonal patch is a 2-connected plane graph in which all inner faces have length 6, 
boundary vertices have degree 2 or 3, and non-boundary vertices have degree 3. These 
graphs are also known as fusenes [18], hexagonal systems [10], poly hexes [17] and 
(6,3)-poly cycles [12] in the literature. These graphs are well-studied in mathematical 
and computational chemistry since they model benzenoid hydrocarbons and graphite 
fragments (see e.g. [18] and the references therein). A central question is that of enu- 
merating hexagonal patches, either of a given size [6], or with a given boundary form. 

A sequence xq, . . . ,JCjt_i of twos and threes is a boundary code of a hexagonal 
patch G if there is a way to label the boundary vertices of G with vq, . . . , v*_i such 
that Vq, . . . , Vfc_i, vq is a boundary cycle of G, and the degree d{yi) = x\ for all ;. Note 
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that cyclic permutations and/or inversions of the sequence can yield different bound- 
ary codes for the same patch, but for the question we study this fact is not important. 
It is well-known and easily observed using Euler's formula that the boundary code of 
a hexagonal patch satifies di — ^3 = 6, where of, is the number of boundary vertices 
with degree ;. We define the parameters daiX) and d^(X) also for sequences X of twos 
and threes, as expected. One may ask whether a hexagonal patch exists that satisfies a 
given boundary code. A result by Guo, Hansen and Zheng [18] shows that this question 
is not as easy as was first expected: in Figure 1 their example is shown which shows 
that different patches may exist with the same boundary code. This can be verified by 
comparing the degree of vi with u\, v% with u%, etc. Our drawing of this graph is taken 
from [8]. (In [18] it is also shown that although multiple solutions may exist, they all 
have the same size.) Therefore the following question should be asked: given a sequence 




Fig. 1. Two different patches with the same boundary code. 



S of twos and threes, how many hexagonal patches exist with boundary code 5? This 
counting problem is called Hexagonal Patch. Guo et al [18] and Graver [17] give con- 
ditions for when solutions are unique, if they exist. Deza et al [10] give an algorithm 
for deciding whether at least one solution exist. The complexity of their algorithm is 
however superexponential. In addition they give a polynomial time algorithm for a very 
restricted case (see Section 3). These results have been generalized to various problem 
variants, mainly by varying the conditions on the face lengths and vertex degrees, see 
e.g. [12,7-9]. However, the question whether the counting problem can be answered 
efficiently remained open. 

In this paper we show that the counting problem Hexagonal Patch can be solved in 
time <9(£ 3 ) where k is the length of the sequence. This is surprising since the number 
of solutions may be exponential in k, as can easily be seen by generalizing the example 
from Figure 1. Therefore, we can only return the number of solutions in polynomial 
time, and not return a list of all corresponding patches. The algorithm can however be 
extended to generate all patches in time • 0(n), where n is the number of returned 
solutions. We remark that it is not hard to generalize our result to the generalizations 
introduced in [8]: A 2-connected plane graph is an (m,k)-patch if all inner faces have 
length k, inner vertices have degree m and boundary vertices have degree at most m. Our 
methods work for instance for (4,4) and (6,3)-patches in addition to (3,6)-patches, but 
for simplicity we restrict to hexagonal patches. 

An additional motivation for this result is the following: fullerene patches generalize 
hexagonal patches by also allowing 5-faces in addition to 6-faces. Such patches model 



fragments of fullerene molecules, and therefore their enumeration is another important 
problem from computational chemistry. Fullerene molecules have at most twelve 5- 
faces. The current result is an essential ingredient for the result we give in a second 
paper [2], where we give a polynomial time algorithm for deciding whether a given 
boundary code belongs to a fullerene patch with at most five 5-faces. 

Our algorithm is based on the following idea: with a few intermediate steps, we 
transform the problem Hexagonal Patch to the problem of counting maximum indepen- 
dent sets in circle graphs. A circle graph G is the intersection graph of chords of a circle 
(detailed definitions are given below). Algorithms are known for the optimization prob- 
lem of finding maximum independent sets in circle graphs [16], but counting problems 
on circle graphs have not been studied to our knowledge. 

In this paper we give a simple dynamic programming algorithm for counting inde- 
pendent sets in circle graphs. In addition this algorithm improves the complexity for the 
optimization problem. Circle graphs can be represented as follows (see Figure 2(a), (b)): 
Every vertex of G is associated with a chord of a circle drawn in the plane, which is a 
straight line segment between two points on the circle, such that two vertices are ad- 
jacent if and only if the two chords overlap (possibly only in a common end). We will 
represent chord diagrams with graphs as follows (see Figure 2(d)). Number the points 
on the circle that are ends of chords with 0, . . . ,k— 1, in order around the circle, and 
view these as vertices. View a chord from i to j as an edge ij. Call the resulting graph 
G' the chord model graph. Note that (maximum) independent sets of the circle graph 
correspond bijectively to (maximum) planar matchings or (M)PMs of G', which are 
(maximum) matchings M that do not contain edges {i,j} and {x,y} with i < x < j < y. 
Hence counting MPMs in G' is polynomially equivalent to counting maximum inde- 
pendent sets in circle graphs. 
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Fig. 2. A circle graph, (simple) chord diagram and (simple) chord model graph. 



Circle graphs are extensively studied and generalize permutation graphs and dis- 
tance hereditary graphs, see e.g. [5]. Recognizing them and constructing a chord rep- 
resentation can be done in polynomial time [4, 15], and the current fastest algorithm 
uses time 0(n 2 ), where n is the number of vertices [22]. A number of problems that 
are NP-hard on general graphs are easy on circle graphs, such as in particular finding 
maximum independent sets [16, 23, 24, 19]. 



The first algorithm for the optimization problem by Gavril [16] has time complexity 
<9(m 3 ), where m is the number of vertices of the circle graph, which is the number 
of edges of the corresponding chord model graph. This was improved to 0{m 2 ) by 
Supowit [23]. Recently this has been improved further by Valiente [24] in the way we 
will explain now. All of these algorithms work with the chord model graph (or chord 
diagram), and as a first step, transform it into a 1 -regular graph as shown in Figure 2(e): 
for a vertex of degree d, d new vertices are introduced, and the d incident edges are 
distributed among these in such a way only one of these edges can appear in a PM 
of G. This does not change the size and number of MPMs. The resulting graph G has 
2m vertices and m edges, and is called the simple chord model graph. We assume the 
vertices are numbered 0, . . . ,2m — 1, in the proper order. Then the length of an edge 
ij G E(G) is \ j — i\. The algorithm from [24] has complexity 0(1), where / is the sum 
of all edge lengths of the simple chord model graph obtained this way. Clearly this is at 
most 0(m 2 ), and in many cases better. However, when dense chord model graphs are 
given on n vertices and m G £2(n 2 ) edges, this algorithm may need £2(n 4 ) steps. Our 
transformation from Hexagonal Patch yields a chord model graph G', which in fact may 
be dense. 

We give a simple algorithm with complexity 0(nm), which not only determines the 
size of a MPM, but also counts the number of MPMs of the chord model graph. This 
improvement in time complexity is possible by working with arbitrary degrees, and 
not using the simple chord model graph, in contrast to all previous algorithms for this 
problem [24]. 

The outline of the paper is as follows. In Section 2 we give definitions, and a precise 
formulation of the problem. In Section 3 we define locally injective homomorphisms 
to the hexagonal lattice (the brickwall) as a way of representing problem instances and 
solutions and reduce the counting problem to a problem on walks in the brickwall. In 
Section 4 we reduce that problem to that of counting proper assignment sets of the walk, 
which is in fact the problem of counting MPMs in chord model graphs. In Section 5 we 
present our algorithm for counting MPMs, and in Section 6 we give a summary of our 
algorithm for Hexagonal Patch. We end in Section 7 with a discussion, where we also 
discuss a similar problem from topology. Statements for which proofs are omitted are 
marked with a star, the proofs appear in the appendix. 

2 Preliminaries 



walk For basic graph theoretic notions not defined here we refer to [11]. A walk of length k 

in a (simple) graph G is a sequence of k + 1 vertices vo, . . . , vj such that v; and v,-+i are 
internal adjacent in G for all i G {0, . . . ,k— 1}. vi, . . . ,\'k-i are the internal vertices and vo, v# the 

end end vertices of the walk. The walk is closed if vo = vj. Throughout this paper we will 

closed in addition assume that v;_i ^ v i+ i for all i G {l,...,k— 1}, and if the walk is closed, 

Vi 7^ V£_i (i.e. we will assume walks do not turn back). If v,- ^ Vj for all i ^ j then 
path the walk is a path. If the walk is closed and v, ^ vj for all distinct i,j G \0,...,k— 1} 

cycle then it is also called a cycle. A cycle of length k is also called a k-cycle. For a walk 

k-cycle W = vo, . . . , V£, W x denotes v x . If W is a closed walk, then W x denotes v x mo( j k- We will 

yj also talk about the vertices and edges of a walk, which are defined as expected. In a 



slight abuse of terminology, the graph consisting of these vertices and edges will also 
be called a walk (or path or cycle if applicable). 

Let H be a hexagonal patch, and B be a boundary cycle of H of length k. Let X = 
xq,... ,Xic-\ be a sequence of twos and threes. We say that the tuple (H,B) is a solution 
for the boundary code X if d(Bj) =Xj for all i G {0, . . . ,k— 1}. Two solutions (H,B) and 
(H' \B') are considered equivalent if there is an isomorphism \\f from H to H' such that equivalent 
yf(Bj) = B'j for all i. Formally, when we ask for the number of different pairs (H,B) that 
satisfy some property, we want to know how many equivalence classes contain a pair 
(H,B) satisfying this property. The counting problem Hexagonal Patch is now defined 
as follows: given a sequence X, how many different solutions (H,B) to X exist? 



3 From Boundary Codes to Walks in the Brickwall 



Fig. 3. The brickwall , 



An (infinite) 3-regular plane graph where every face has length 6 is called a brick- 
wall. It can be shown that the facial cycles are the only 6-cycles of a brickwall, and that 
all brickwalls are isomorphic. 

We will use SB to denote the brickwall as drawn in Figure 3. Edges that are hori- 
zontal (vertical) in this drawing are called the horizontal (vertical) edges of 3§. Paths 
consisting of horizontal edges are called horizontal paths. Two vertices joined by a 
horizontal path are said to have the same height. 

The reason that we study brickwalls is because the following mapping of hexagonal 
patches into them is very useful. Let H be a hexagonal patch. A locally injective homo- 
mophism (LIH) of H into S§ is a mapping of the vertices of H to vertices of 33, such 
that adjacent vertices are mapped to adjacent vertices in S3, and such that all neighbors 
of any vertex in H are mapped to different vertices in S3. Since the shortest cycles in SB 
are of length 6, a LIH into 23 maps 6-cycles to 6-cycles. Since the faces of 38 are the 
only 6-cycles in 33, we see that a LIH of H into 33 also maps inner faces to faces. 

Loosely speaking, the idea behind these mappings is as follows. Let H be a hexag- 
onal patch of which we fix a boundary cycle B. When we map H with a LIH into 3S, 
then the boundary B is mapped to some walk W in 33. But now it can be shown that 
this walk W is only determined by the choice of the initial vertices and the boundary 
code of H. Hence instead of asking how many hexagonal patches exist with a certain 
boundary code, we may ask how many patches exist that can be mapped properly to 
the brickwall, such that the boundary coincides with the walk that is deduced from the 
boundary code. Below we will go into more detail. 
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The technique of mapping patches to brickwalls is not new, and is actually consid- 
ered folklore to some extent [10]. For instance, Deza et al [10] observe that Hexagonal 
Patch can be solved in polynomial time if the LIH is bijective, and Graver [17] shows 
that the problem Hexagonal Patch can only have multiple solutions if there is a brick- 
wall vertex that has at least three preimages in such a LIH. We will however study these 
mappings more in more detail than has been done before, and develop new concepts, 
and prove new statements which we feel are of independent interest. 

Let W be a walk in a 3-regular plane graph G. We say W makes a right (left) turn at 
i when edge WJWi_i immediately follows edge W,W, + i in the clockwise (anticlockwise) 
order around W,. Note that since we assume that walks do not turn back and G is 3- 
regular, W makes either a left or a right turn at every i. 

Walk construction: Using a given sequence xq, . . . of twos and threes, we con- 
struct a walk W = Vo,. . . in 3§ as follows. For Vo and vi, choose two (arbitrary) 
adjacent vertices. For i > 1, choose v,-+i such that W makes a left turn at ; if x, = 3, and 
makes a right turn at i if Xj = 2. 

Let W be a closed walk in 33 of length k, H be a hexagonal patch, a LIH from 
H to 33 and B a boundary walk of H of length k. Then the tuple (H,(j),B) is said 
to be a solution for W when 0(B,-) = W, for all i. Two solutions S = (H,<j>,B) and 
S' = (H' , <j)',B') are considered to be equivalent if and only if there is an isomorphism 
y from H to H 1 such that yf(Bj) = B\ for all i. We say that y/ is an (or demonstrates the) 
equivalence between S and S'. The LIH allows us to use the terminology defined for 
3§ for the graph H as well; we will for instance call edges of H horizontal or vertical if 
their images under are horizontal or vertical, respectively. 

Let the boundary B of a hexagonal patch H be mapped to the closed walk W in 38 
by the LIH 0. This is a clockwise solution if and only if for every i, d{Bi) = 2 if W 
makes a right turn at i, and d(Bj) = 3 if W makes a left turn at i. It is anticlockwise 
when these conditions are reversed. Let Right(W) and Left(W) denote the number 
of indices i S {0, . . . , k — 1 } such that W makes a right turn or left turn at i, respectively. 
The turning number of W is t(W)= (RIGHT (W) - Left(W))/6. Using the fact that 
for a solution (H,(j>,B), (j) maps faces of H to faces of 38, it can be shown that every 
solution is either clockwise or anticlockwise. Since a hexagonal patch has di — d?, =6 
(di is the number of degree i vertices on the boundary), Lemma 1 then follows. Variants 
of Lemma 2 have been proved in [8, 17]. 

Lemma 1 (*) Let W be a closed walk in 3§. If t (W) = 1, then every solution to W is 
clockwise. Ift(W) = — 1 then every solution to W is anticlockwise. Ift(W) $ { — 1, 1}, 
then no solution exists. 

Lemma 2 Let (H,B) be a solution to a boundary code X and let W be a walk in 3§ 
that is constructed using X. Then there exists a unique LIH <j) such that (H,(j),B) is a 
clockwise solution to W. 

Because of Lemma 2, we may rephrase the problem Hexagonal Patch in terms of 
solutions (H,(p,B) to a closed walk W in the brickwall. 

Theorem 3 The number of different (hexagonal) solutions for a boundary code X with 
^2p0 — dj,(X) = 6 is the same as the number of different clockwise solutions for the 
walk W in S3 that is constructed using X. 



Proof: For any solution (H,B) for X, a unique LIH exists such that (H,<j>,B) is a 
clockwise solution to W (Lemma 2). For any clockwise solution (H,<j>,B) to W, the 
characterization of clockwise solutions and the construction of W shows that (H,B) is 
a solution to X. (Since d2(X) — dj(X) = 6 and t(W) = 1 by Lemma 1, W turns at as 
prescribed by xo.) Note that the definitions of equivalence for pairs (H,B) and triples 
(H,(j),B) coincide and, in particular, do not depend on (j). □ 



4 From Walks in the Brickwall to Assignment Sets 



Throughout Section 4, W denotes a closed walk in with length k. We first sketch the 
main idea of this section. If we consider a solution (H,<j>,B) to W, then we mentioned 
above that this defines which edges of H are horizontal and vertical. Now if we start at 
a boundary vertex Bj of H that is incident with a horizontal interior edge of H, then we 
can continue following this horizontal path of H until we end in a different boundary 
vertex Bj. We will say that this solution assigns i to j. If we only know all assignments 
defined by the solution this way, we can reconstruct the unique solution. We will deduce 
properties of such sets of assignments such that there is a solution if and only if these 
properties are satisfied. The purpose is to show that we may focus on counting such 
assignment sets instead of solutions to the walk. 

For all i, j where W, and Wj lie on the same height, H, y denotes the horizontal walk 
in 3& from W, to Wj. Consider an index / G {0, 1} and the vertex Wj. Let u be 

the neighbor of Wj in 38 not equal to Wi-\ or W;+i. If u has the same height as Wj and 
W makes a left turn at i, then index i is called a PA-index. In Figure 4(b) an example 
is shown, where vertices corresponding to PA-indices are encircled, and their indices 
are shown. Note that if W has a clockwise solution (H,(j),B), then the PA-indices are 
precisely those indices i such that Bj has degree 3 and the interior edge incident with Bj 
is horizontal (see Figure 4(a)). 

A possible assignment (PA) is a pair {/, j} of PA-indices with Wj ^ Wj such that Wj 
and Wj have the same height and Hjj does not contain any of Wj-i, Wj+i,W/-i, Wj+i 
(note that Hjj has non-zero length). For instance, in Figure 4(b) some PAs are { 1 , 28}, 
{ 1 , 14} and {21 , 32}, but { 1 , 24} is not. 
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(a) Assignment paths in the patch H (b) PA-indices of the walk W 

Fig. 4. Assignment paths of a solution and PA-indices of a walk. 



Let (H,<j),B) be a clockwise solution to a closed walk W in S3. An assignment path 
P is a horizontal path in H from B, to Bj where ; ^ j, and all edges and internal vertices 
of P are interior edges and vertices of H. In Figure 4(a) the assignment paths of the 
given solution are shown in bold. 

Proposition 4 (*) If a clockwise solution (H,<j),B) to W contains an assignment path 
from Bi to Bj, then {i, j} is a PA ofW. 

This motivates the following definition. A clockwise solution S = (H,<j>,B) to a 
walk W assigns i to j if there is an assignment path from Bj to Bj. For each clockwise 
solution 5, we define the set si (5):= {{i,j} '■ is a PA of W and S assigns i to j}. 
This is the assignment set defined by the solution S. 

Lemma 5 (*) Let W denote a closed walk in 3§ and let 5,5' be clockwise solutions of 
W. IfS and S' are equivalent, then si (5) = si (5'). 

Now we will deduce the properties of a set si (5). Proposition 6 shows that assign- 
ment paths do not share vertices. Combining this with planarity yields Proposition 7. 

Proposition 6 (*) Let (H,<j),B) be a clockwise solution to W. Every interior vertex of 
H and every vertex Bi, where i is a PA-index, lies on a unique assignment path. 

Proposition 7 (*) Let S be a solution to a closed walk W that assigns i to j. For any 
x,y with x < i <y < j or i < x < j < y, S does not assign x to y. 

These two propositions give us properties a set of the form si (5) for a clockwise 
solution 5 necessarily has to have. Given W, a set A of possible assignments of W is a 
perfect matching on the set of PA-indices if for every PA-index ; of W there is exactly 
one pair {i,j} E A. A is non-crossing if there do not exist assignments {i,j},{x,y} G A 
such that i < x < j < y. An assignment set for W is a set of possible assignments of 
W. It is a proper assignment set if it is a non-crossing, perfect matching on the set 
of PA-indices of W. Combining Proposition 4, Proposition 6 and Proposition 7 yields 
Lemma 8. Lemma 9 states more or less the reverse; the long proof appears in Ap- 
pendix B. 

Lemma 8 IfS = (H, (j>,B) is a clockwise solution ofW then si(S) is a proper assign- 
ment set for W. 

Lemma 9 (*) Let W denote a closed walk in SS with t(W) = 1, and let A be a proper 
assignment set ofW. Then there exists a clockwise solution S ofW with si{S) = A. 

It remains to establish the converse of Lemma 5. Suppose we have two solutions 
S = (H,(j),B) and S' = {H',(j)',B') with si(S) = si(S'). Every vertex of H and H' lies 
on the boundary or on an assignment path (Proposition 6). Therefore we can use the 
boundary and the assignment paths to define a bijection y/ : V(H) — ► V(H'). When 
doing this appropriately, it can be shown that y an equivalence. 

Lemma 10 (*) Let W be a closed walk in 38, and let 5 and S 1 denote clockwise solu- 
tions ofW. If si (5) = si (5'), then S and S' are equivalent. 



Theorem 11 Let W be a walk in 3$ with t(W) = 1. The number of equivalence classes 
of solutions to W is the same as the number of different proper assignment sets for W. 



Proof: The above lemmas show that S i— > s/ (S) gives a bijection from the set of equiv- 
alence classes of clockwise solutions of W to the set of proper assignment sets for W, 
since the following properties are satisfied: (1) s/ is well-defined: Let S\ and S2 denote 
clockwise solutions of W. If Si and S2 are equivalent, then s/(S\) = s/ (S2) (Lemma 5). 
(2) The range of s/ is correct: For any clockwise solution S of W the set sf (S) is a 
proper assignment set for W (Lemma 8). (3) s/ is injective: Let Si and S2 denote clock- 
wise solutions of W. If sf(S\) = si [Si), then Si and S2 are equivalent (Lemma 10). 
(4) si is surjective: For any proper assignment set A for W, there exists a clockwise 
solution S of W with sf(S)=A (Lemma 9). □ 

It follows that for solving the Hexagonal Patch problem, we may focus on counting 
proper assignment sets for the walk W (assuming t(W) = 1). 

5 Counting Maximum Planar Matchings 

In this section we will observe that the remaining algorithmic problem is that of count- 
ing independent sets in circle graphs, and present a fast algorithm for this problem. We 
use the closed walk W in 88 to construct a graph G w with vertex set V = {0, . . . , n — 1}, 
where n is the number of PA-indices of W. Let po,... ,p„-i be all PA-indices of W, 
numbered according to their order in W. Then the edge set of G w will be E = {ij 
{pi,Pj} is a PA of W}. The following lemma is now easily observed. 

Lemma 12 Let G w be the graph as constructed above from the walk W. IfG w has no 
perfect planar matching, then W has no proper assignment sets. Otherwise the number 
of proper assignment sets for W is equal to the number of MPMs in G w . 

Now we will present an algorithm for counting MPMs of a graph G with V(G) = 
{0, ...,« — 1}. As mentioned in the introduction, this is equivalent to counting maxi- 
mum independent sets in a circle graph H, where G is the chord model graph of H. We 
will present this algorithm for the general case where G has edge weights: wij denotes 
the edge weight of ij, and a PM M is maximum if Y*eeM w e is maximum. 

For i,j £ V(G) with i < j, let G,;= G [{i, ... ,j}]. If i > j, then Gu is the empty G/j 
graph. For i,j £ V(G), let Sy denote the size of a MPM in Gy. In particular, So.n-i Si j 
is the size of a MPM in G. We now give a subroutine S(i,j) for calculating S,j, which 
considers the sizes of various PMs for G,j, and returns the size of the largest PM. 

A subroutine S(i,j) for calculating Su: 

(1) m:=0 

(2) If r < thenm := S i+ ij 

(3) For v e N(i) with i + 1 < v < j - 1: 

(4) m := max{/n,S !iV + 5v+ij} 

(5) If j £ N(i) and ; > i then 

(6) m: = max{m, wij + S,+ 1 j- 1 } 

(7) Return m 



Lemma 13 (*) Let G be a graph with V(G) = {0,...,n — 1} and i,j G V(G). If the 
values S X j are known for all x,y with —I <y—x< j — i, then the subroutine S(i,j) 
computes Sij in time 0{d(i)). 

Let Nij denote the number of MPMs in G; ,-. Below is a similar subroutine N(i,j) 
for calculating Nu, which considers various PMs for G,- ,-, checks whether they are 
maximum by comparing the size with Sij, and keeps track of the number of MPMs 
using the variable N. 

A subroutine N(i,j) for calculating N,j: 

(1) If j </ then Return 1, exit. 

(2) N:=0 

(3) If Stj = S MJ then N :=N+N MJ 

(4) For v G N(i) with i + 1 < v < j - 1: 

(5) If S u = S iA , + S v+Lj then N :=N + N Lv x N v+lJ 

(6) If ; G N(i) and 5/j = w y + then A? := N+N i+1 j-i 

(7) Return^ 

Lemma 14 (*) Lef G be a graph with V(G) — {0, . . . ,n - 1} anc/ i,j G V(G). T/'f/je 
values S x ,y and N x> y are known for all x,y with — 1 < y — x < j — i, and Sij is known, 
then the subroutine N(i,j) computes Nij in time 0(d(i)). 

Theorem 15 Let G be a graph with V(G) = {0, ...,« — 1} on m edges. The size and 
number of MPMs of G can be computed in time 0(nm). 

Proof: For d = — 1 to n — 1, we consider all i,j G {0,...,n — 1} with j — i = d, and 
calculate Sij and Nij using the above subroutines. This way, for every value of d, every 
vertex of G is considered at most once in the role of i. For this choice of i, calculating 
Sij and Nij takes time 0(d(i)) (Lemma 13, Lemma 14). Hence for one value of d this 
procedure takes time 0(L lG y(G) d (/)) = 0{m). □ 

We remark that Valiente's algorithm [24] for simple (1 -regular) chord model graphs 
can also be extended by using Subroutine N(i,j) to calculate Nu in constant time, 
immediately any time after a value 5,- ; is calculated. This then yields time complexity 
0(1) and space complexity 0(n). In some cases it may be better to transform to a simple 
chord model graph and use this algorithm. 

6 Summary of the Algorithm 

We now summarize how counting the number of hexagonal patches that satisfy a given 
boundary code X of length k can be done in time 0(k 3 ). W.l.o.g. d 2 (X)-d 3 (X) = 6. 
First use X to construct a walk W in 88 of length k, as shown in Section 3. Theorem 3 
shows that we may now focus on counting clockwise solutions to W. If W is not closed 
it clearly has no solution. Since d2(X) — d^(X) = 6 we may now assume t(W) = 1. 
Then Theorem 1 1 shows we may focus on counting proper assignment sets for W. Now 
construct G w as shown in Section 5. G w has n vertices where n < k is the number of 
PA-indices of W (and 0(n 2 ) edges). By Lemma 12, the number of proper assignment 
sets for W is equal to the number of MPMs of G w , provided that G w has a perfect PM. 
This number and property can be determined in time <9(n 3 ) G 0(k 3 ) (Theorem 15). 



7 Discussion 



Our first question is whether the complexity of 0{k i ) can be improved. Secondly, con- 
sidering the motivation from benzenoid hydrocarbons, it is interesting to study whether 
a patch exists that has a 'reasonably fiat' embedding in R 3 using regular hexagons. 
More precisely, this is the brickwall walk problem, but requires in addition giving a 
consistent linear order ('depth') for all vertices mapped to the same vertex of 38. It may 
also be interesting to study generalizations such as to surfaces of higher genus. 

After we presented an early version of this work [3], Jack Graver pointed us to a 
similar well-studied problem in topology. Let S l denote the unit circle and D 2 the unit 
disk in R 2 . An immersion is a continuous function / : A — > B such that for every x in 
A there is a neighborhood N of x such that is a homeomorphism. (A curve when 
A = S , B = R 2 .) An immersion c : S l — » M 2 of the circle into the plane is normal 
if c has only finitely many double-points and c crosses itself at each of these. Two 
immersions d, d 1 are equivalent if there exists a homeomorphism <j) : R 2 — > IR 2 such 
that d o = d' . Now the Immersion Extension problem is this: given an immersion 
c : S l — > R 2 , how many immersions d : D 2 — > R 2 exist that extend c? Note that this 
problem is not combinatorial, therefore it makes no sense to study its computational 
complexity. One can turn it into a combinatorial problem by restricting the input to 
piecewise linear (PL) curves c : S l — > R 2 . 

When viewing the walk constructed in Section 3 as a curve, there are obvious sim- 
ilarities between the Hexagonal Patch problem and the Immersion Extension problem. 
However, to our knowledge it is an open problem to prove that these problems are in 
fact equivalent. The ideas introduced here may be helpful for giving such a proof. Es- 
tablishing this would provide insight to both problems, since the Immersion Extension 
problem is well-studied - at least on normal curves - see e.g. [1, 13, 21]. Interestingly, 
Blank [1, 14] reduces the Immersion Extension problem problem to a combinatorial 
problem that is essentially the same as counting MPMs in simple chord model graphs. 
He does not address the complexity of this problem. Shor and Van Wyk [21] were the 
first to study the complexity of the combinatorial Immersion Extension problem on 
normal curves. They give an (9(« 3 log«) algorithm where « is the number of pieces of 
the PL curve c. Assuming the equivalence of the Immersion Extension problem and 
the Hexagonal Patch problem, this would give an alternative algorithm for Hexagonal 
Patch; note that there are methods for transforming general PL curves to equivalent nor- 
mal PL curves [20]. Since our algorithm does not need such a step, it is not only faster 
but also much easier to implement (see also [20]). However, the question of equivalence 
of these problems is still interesting because many generalizations of the Immersion Ex- 
tension problem have been studied [13]. Finally, we believe that in fact our method can 
be adapted to give a simple and fast for the combinatorial Immersion Extension problem 
that does not require the assumption that the given curve is normal, but that is beyond 
the scope of this paper. 
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A Proofs of Section 3 and Section 4 



Proof of Lemma 1 : Let be a LIH from a hexagonal patch H with boundary cycle 
B of length k to SB, and let W = <j)(B ), 0(Bi), . . . , 0(5jt-i), 0(J5 O ). We first show that 
(H,<j>,B) is either a clockwise solution to the walk W or an anticlockwise solution to 
W. 

An index i is called locally clockwise if either d(Bj) = 2 and W makes a right turn 
at i, or c/ (B;) = 3 and W makes a left turn at i. 

We show that if some i is locally clockwise, then every index is locally clockwise. 
Suppose this is not true, so then there is an i that is locally clockwise such that i + 1 
is not. Assume first that d(Bj) = 2 and c/(B, + i) = 2. Then W makes a right turn at i, 
but a left turn at i + 1. Therefore W/_ i and W/+2 do not lie at a common facial cycle of 
SB. Since B, and B 1+ i both have degree 2 in H, all of the vertices B,_i , . . . ,B,+2 lie at a 
common inner face of H. This is a contradiction since <j> maps faces of H to faces of SB. 

In the case where d(B,) = d(Bj + i) = 3, we consider the neighbor v, of B,- that is 
not equal to B,_i or B, + i, and the neighbor v, + i of B, + i that is not equal to B, or B, + 2. 
These again lie at a common face of i/, but if W makes a left turn at i and a right 
turn at /+ 1, are mapped to two vertices that do not lie at a common face, which again 
yields a contradiction. The two other cases are analogous. We conclude that if a solution 
contains a locally clockwise vertex, it is clockwise. 

Now we relate this to the turning number. Let c/, denote the number of vertices of 
degree i on the boundary of H. If (H,(j>,B) is a clockwise solution then d^ = Right(W) 
and c/3 = Left(W). We know that di — d^ =6 since H contains no 5-faces. Hence 
t(W) = 1. Similarly, if an anticlockwise solution exists then t(W) = — 1 follows, which 
proves the statement. □ 

Proof of Proposition 4: We show that if a clockwise solution (H,<j>,B) toW assigns 
Bi to Bj, then {/, j} is a PA of W. Let P be an assignment path from 5,- to Bj. We have 
that Bj and Bj have the same height, since P is horizontal. B, ^ B ; holds since H is 
2-connected. P is then mapped to a non-zero length path in 8§ (it does not turn back, 
since is a LIH), so Wi 7^ Wj follows. All edges of P are interior edges of H, so B; and 
Bj have degree 3, and therefore W makes a left turn at i and j. It follows that ; and j are 
PA-indices. Since <p is a LIH and P contains no boundary edges, H, j does not contain 
any aSWi-uWi+uWj-uWj+t. □ 

Proof of Lemma 5: We show that if two clockwise solutions S and S 1 of W are 
equivalent, then sf(S) = af{S'). Let S = (H,<j),B) and S' = (H' ,<j>' ,B'), and let y : 
V(H) — > V(H') demonstrate their equivalence. Note that <p' o y/ and (j) are both LIHs 
from H into SB that map B; to Wj. As the LIH with this property is uniquely determined 
by Lemma 2, we conclude that <j>' o = 0. In particular, any edge mv is horizontal in 
// if and only if y/(u)y(v) is horizontal H' . Clearly an analogous statement holds for 
vertices being interior. Therefore y maps assignment paths to assignment paths. Since 
y{Bi) = B\ for all i it follows that srf(S) = #/(S'). □ 

Proof of Proposition 6: Let (H,<j>,B) be a clockwise solution to W. We show that 
every interior vertex of H and every vertex B,, where i is a PA-index, lies on a unique 
assignment path. Let M be the set of horizontal non-boundary edges of H, and let H' = 



(V(H),M). Since is a LIH, H' has maximum degree at most 2. //' contains no cycles, 
because these would have to be mapped to cycles of 3S but S3 contains no cycles with 
only horizontal edges. Hence H 1 is a set of paths and isolated vertices. 

It can be seen that vertices with degree 2 in H' are interior vertices of H, and that 
vertices with degree 1 in H' are equal to B, for some PA-index i. Hence the path com- 
ponents of H' (paths of non-zero length) are all assignment paths. Since all assignment 
paths in H are also part of H', we see that there is a one-to-one correspondence between 
assignment paths in H and non-trivial components of H' . We also see that every inte- 
rior vertex of H and every vertex B, where i is a PA-index lies on one such path. The 
statement follows. □ 

Proof of Proposition 7: We show that if a solution (H,(j) 7 B) assigns i to j and i < 
x < j < y, then it does not assign x to y. Suppose H contains an assignment path P from 
Bi to Bj, and an assignment path Q from B x to B y . By Proposition 6, P and Q have 
no vertices in common. But since the (distinct) end vertices of the paths appear in the 
order Bi, B x , Bj and B y along a boundary cycle of the plane graph H, this is impossible. 
(Formally, to obtain a contradiction, we may use P, Q and the boundary cycle of H to 
exhibit a subdivision of K4 that is embedded with all vertices on the boundary, which 
then would yield a planar embedding of ^5.) □ 

For the proof below and later proofs in the appendix, it is important to distinguish 
horizontal left between two different kinds of horizontal edges of SB: horizontal left {horizontal right) 

horizontal right edges are edges that follow a vertical edge after turning left (right). Note that this par- 

titions the edges of SB into vertical edges, horizontal left edges, and horizontal right 
edges, and that every face contains two of each. The same holds for faces in a solution 
(H,(j),B); Recall that if patch H is mapped by a LIH to SB, this allows us to define 
vertical and horizontal (left / right) edges in H. Similarly, we will talk about vertices of 
H that lie to the left / below etc. other vertices. This is also defined by (j) and the chosen 
drawing of &8. 

Proof of Lemma 10: We show that if s/{S) = si (S') for two solutions S and S', then 
these solutions are equivalent. Let S = (H,<j>,B) and S' = (//', <j>',B'). We construct the 
isomophism yr from H to H 1 that will demonstrate the equivalence as follows. For all 2, 
y/(Bi) = B\. This defines l// for boundary vertices. Every non-boundary vertex lies on a 
unique assignment path (Proposition 6). Suppose such a vertex v lies on an assignment 
path P from Bj to Bj. Then {/, j} £ srf (S) and thus {;, j} £ stf (S'). The assignment path 
P' from B\ to B'j in H' is also mapped by (j)' to Htj and therefore has the same length 
as P (since (j) and <p ! are LIHs). Now if v is the x-th vertex on P, y/ will map v to the 
x-th vertex of P' . This defines Since every vertex of H lies on the boundary or on 
an assignment path, the function y/ is defined for every vertex of H, and since the same 
holds for H' , \j/ is a bijection. By definition y/ maps boundary vertices to the correct 
boundary vertices, so to demonstrate that yr is an equivalence between S and S', it only 
remains to show that it is an isomorphism. 

We only show that edges of H are mapped to edges of H' by yt. By symmetry a 
similar statement then follows for yr~ 1 , which proves that yr is an isomorphism. Clearly 
y> maps boundary edges of H to boundary edges of H'. Observe that every horizontal 
non-boundary edge of H lies on an assignment path. Therefore yt also maps horizontal 



edges of H to edges of H'. What remains are vertical edges of H that do not lie on the 
boundary. Note that y/ maps horizontal left (right) edges of H to edges of H' of the 
same type, and by observing the same for y/^ 1 , it also follows that if a vertical edge is 
mapped to an edge, it is mapped to a vertical edge again. 

Suppose there exists a (vertical, interior) edge of H that is not mapped to an edge 
of H 1 by y/. Let e = u\\'\ G E(H) be such an edge such that all edges that lie to the left 
of it are mapped to edges of H' . Suppose u\ lies below vi. e is incident with two inner 
faces of H, so we may choose F = u\,U2,U3,V3,V2,V\,u-[ to be the inner face of H on 
the left side of e. Note that all other edges of F are mapped to edges of H'; four edges 
are horizontal, and the other vertical edge is mapped by our choice of e. 

Let u\, u' 2 , w 3 , v 3 , v' 2 and v[ respectively be the images under 1// of the vertices of F. 
The edges v\v' 2 and v 2 v 3 lie on a common face of H' (since v 2 has degree at most 3). 
We show that they lie on a common inner face of H' . If not, then both v \ v 2 and v 2 v 3 are 
boundary edges. Then the corresponding edges V1V2 and V2V3 of H are boundary edges 
too, and thus these two edges share both an inner face and the outer face. Since H is 
2-connected, it follows that d (V2) = 2. Because both H and H' are clockwise solutions, 
d(v 2 ) = 2. Hence the two edges in H' also share two faces, and thus one inner face. 

Let F' = v'j, v' 2 , v' 3 ,x,y,z, v\ be the inner face of H' on which these two edges lie. 
Since \j/ maps to edges of the same type, v\v' 2 and v 2 v 3 are horizontal left and hori- 
zontal right edges respectively. Since 0' maps inner faces of H' to inner faces of 38, 
v'^x is a vertical edge. v' 3 is incident with at most one vertical edge (0' is a LIH), 
so we may conclude that x — m 3 (here we use the fact that vertical edges are not 
mapped to horizontal edges, so v 3 zt 3 is vertical). Continuing this reasoning shows that 
F' = v\,v' 2 ,v' 3 ,u' 3 ,u 2 ,u' i ,v' i . Hence u\v\ G E{H'), a contradiction with the choice of 
mi vi . We conclude that y/ is an isomorphism, which concludes the proof. □ 



B The proof of Lemma 9 



Before we can prove Lemma 9 we need to introduce some new terminology and lem- 
mas. For a closed walk W = Vo, • • • ,V£, the subwalk of W from i to j is the walk 
v/, v,-+i, . . . ,vj of length j — i mod k. If j < i then, more precisely, this is the walk 
v,-, v ( +i, . . . ,V£_i,Vo, . . . ,v,-. The subwalk of W from i to j will be denoted by Wy. Wjj 
If W = v'o, . . . ,v# is a closed walk, then for any i 6 {0,...,k— 1}, the walk W' = 
Vj, v,-+i, . . . , V£_i, vo, vi, . . . ,Vj-i, v, is called a rotation of W. We will write W' ps W to 
express that W' is a rotation of W. For a pair of walks W = vo, . . . ,v# and W' = uq, . . . ,«/ 
with Vfc = «o an d v k-\ 7^ M i> W o W' denotes the concatenation of W and W*', which 
is vo, . . . ,V£,wi, . . . ,Ui. In sequences, the notation {a) b means that b copies of a are 
inserted in the sequence at this point. For instance, 1,2, (3) 3 , 4 denotes the sequence 
1,2,3,3,3,4. 

Let W be a closed walk in 38. For indices i, let « (?) be the first PA-index after i (not 
equal to i). So with respect to the walk W shown in Figure 4(b), «(1) = 4, «(32) = 0, 
etc. The subwalks of the form W; for any PA-index i are called the pieces of W. pieces 



rotation 
WoW' 

concatenation 
{a? 

n(i) 



Elementary Cycles Let A be a proper assignment set of a closed walk W in 38. We 
will now define how such a tuple W,A gives elementary cycles, which correspond to 



closed walks in SB. Informally, for any index i the unique elementary cycle that con- 
tains W/W/+1 can be found as follows. This is illustrated in Figure 5, where the proper 
assignment given by the solution shown in Figure 4(a) is used. Arcs are shown to in- 
dicate the direction of the elementary cycles, and the head of the arc indicates the first 
vertex. 
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(a) Walk W and proper assignment A 



Fig. 5. Elementary cycles in 3§ gi 



(b)The elementary cycles of 
by a proper assignment A of W. 



Start at W;, and continue along W (in the direction of increasing indices) until the 
first PA-index n(i) is found. Let {n(i), i\} 6 A be the (unique) assignment in A that 
contains «(/). Then continue along the horizontal path #„(,), ;i in SB. At Wi lt follow 
W again until the next PA-index n(ii), and then follow H n (i{),ii where {n(ii),i2} G 
A, etc. Continue with this procedure, taking alternatingly pieces of W and horizontal 
paths that correspond to assignments in A until we arrive again at Wj. Observe that this 
procedure ends and that we do actually arrive at Wj again this way. We consider the 
choice of starting vertex to be irrelevant, hence in the following formal definition we fix 
elementary cycle a canonical rotation. An elementary cycle of W, A is a walk in SB of the form 

C = W kAio] o H n[k)M oW hMh] o . . .oH n(il)tk , 

where {n(ij),ij+i} E A for all j, and {«(//), io} £ A. In addition, we require that zq < ij 
for all j G {1,...,/}. This last condition fixes the canonical rotation. Note that ele- 
mentary cycles that start at different indices of W may still yield the same walk in SS. 
Therefore we also consider a sequence / of numbers that give the corresponding in- 
dices of W. We insert —1 in this sequence for the vertices that correspond to vertices of 
horizontal paths instead of parts of W. For example, for the above choice of C, 

/ = /o,.. .,«('o) . (-i)* -1 i h,---Mh) , , , h), 

where x and y are the length of ff„(j ) ^ and H n / t \ iQ respectively. (Obviously, if iq > 
n(io), then io, . . . ,n(z'o) should be read as j'o , i'o + 1,0,1, ... ,«(*'o)> where k is 

the length of W, etc.) 

Formally, an elementary cycle is now a pair C,I of a closed walk C in S§ and se- 
quence of numbers / that are of the form explained above. This formal definition is 
needed to clearly define what the number of elementary cycles of W,A is: elementary 



cycles are still considered different even if they yield the same walk, but admitting 
different rotations is irrelevant. However, below we will often informally denote ele- 
mentary cycles just by C; the index sequence / is clear from how we denote C. Note 
that for every i there is a unique elementary cycle C,I such that / contains i and i + 1 
consecutively. In a slight abuse of notation, from now on we will often simply express 
this statement as follows: there is a unique elementary cycle C that contains the walk 
edge WjWj + i. In the special case where A = 0, there is only one elementary cycle C,/, 
which has C = W. 

A walk in 8$ in which every vertical edge is followed by a horizontal right edge and 
preceded by a horizontal left edge is called a right-turn walk. right-turn walk 

Proposition 16 Let A be a proper assignment set for a walk W in SB. Then every ele- 
mentary cycle C ofW,A is a right-turn walk. 

Proof: Any vertical edge of C must come from a piece of W. So let W/_iW,- be this 
vertical edge. If W turns left at i, then i is a PA-index and C turns right. If W turns right 
at i, then i is not a PA-index and C turns right as well. Hence in both cases, a vertical 
edge in C is followed by a horizontal right edge. Similarly, C turns right at i — 1 in both 
the case that it is a PA-index and the case that it is not, hence vertical edges in C are 
preceded by horizontal left edges. □ 

For the following proof, we use the following vertex labelling for SB. See also Fig- 
ure 6, which illustrates the next Proposition. 

V(,3B) = {b u :iJeZ} 

E(SB) = {bijb i+lJ : i,j 6 Z} U { /'/. ;/■>< i . , . i : ij £ Z,i odd} 



bi.i 



''4.1 



W, 
Wo 



''5.1 
''4.0 



''7.1 
''(,.(1 



bs.u I *5.0 I ^7.0 

Fig. 6. A right-turn walk W in 3§ with t(W) = 2. 



Proposition 17 Let W be a closed right-turn walk in 33. Then t(W) > 1. IftiW) = 1, 
then W has a clockwise solution. 

Proof: For a closed walk W in SB of length k and any / £ {0, . . . ,k} we define t(W,l) 
to be the number of indices ; with 1 < i < I such that W makes a right turn at i minus the 
number of those indices where W makes a left turn. So t (W, 0) = and t (W, k) = 6t(W), 
and for every I, t (W, I + 1) = t (W, I) ± 1 . 

Every closed walk in SB contains a vertical edge (since walks do not turn back), so 
w.l.o.g. assume WqWi = b\ t 0^2.1- Since W is a right-turn walk, after this a horizontal 



right edge follows, which is part of an alternating sequence of horizontal right and left 
edges. This sequence continues until a horizontal left edge is followed by a vertical 
edge, and after that an alternating sequence of horizontal right and left edges again 
follows (note that the walk cannot close before this point). So for some i > 1 and j > 1 
we have 

Wo,2i+2y+2 = ^1,0)^2,1)^3,1) • ■ • )&2H-2,l>&2j+l,0>&2(,0, ■ ■ ■ ,&2i-2,/+l,0- 

(In Figure 6, / = 2 and j — 1. Note that in general j > i is also possible.) Choose j 
maximum, so either the walk W closes at this point (k = 2i + 2j + 2), or W continues 
with another vertical edge. In either case, the sequence t(W,0),t (W, 1), . . . ,t (W,2i + 
2j + 2) is then 0, (1,2)', 3, (4,5) ; ,6. Continuing this reasoning shows that t(W,l) can 
never decrease below 6 when I > i + j + 2, so we conclude t(W) = t(W,k)/6 > 1. In 
addition, if t(W) = 1, then the walk cannot contain another vertical edge, so 

W = ^1,0,^2,1)^3,1,- ■ ■ ,£>2/+2, 1,^2/4- 1,0, ^2/,0, ■ ■ -,&1,0- 

It is easily seen that in this case the subgraph of induced by the vertices of W is a 
clockwise solution to W. □ 
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Fig. 7. Splitting the walk IV and assignment set A into two. 



Splitting a walk In the following lemmas we use the following notations. Let A again 
be a proper assignment set for the closed walk W in S3. We assume A contains a PA 
{Q,y} such that n(y) = (i.e. W y o is a piece of W). We now consider the two new closed 
walks W'— Wo, y oHyfl and W"= W v ,o °H 0>y in ^. Note that W" is in fact an elementary 
cycle of W,A, and that for all ; £ {0, . . . W/ = W,-. Figure 7 illustrates this for the 
walk and assignment from Figure 5 (note that y = 32). 

We remark that although we consider various walks in S3 in this section, the nota- 
tions Hjj and «(/) are always defined with respect to W. Let A'= A — {0,y}. (We denote 
A\{x} and A U {x} as A — x and A +x respectively.) Observe that i is a PA-index of W' 
if and only if it is a PA-index of W and ; ^ 0,y. Thus it is easily seen that: 



Proposition 18 A' is a proper assignment set for ff'. 



The above construction of smaller walks and corresponding assignment sets from a 
given pair ff, A allows various induction proofs, justified by the next lemma. 

Lemma 19 Let W, A, ff ' and A' be as defined above. Then 

1. For every elementary cycle C',I' ofW' ,A', there is an elementary cycle C,I ofW,A 
with C « C. 

2. W,A has one more elementary cycle than W' ' ,A' '. 

Proof: We use the notation n' and H[ j for the walk W', which are similar to the notations 
n and //,■ j for ff , so n'(() denotes the next PA-index of ff ' after i, and H-j denotes the 
horizontal path between W( and Wj. Let C',I' be an elementary cycle of W',A' with 

C' = ff' a- \ oH 1 a- n ■ off. ,r. <. o ... oH 1 a - \ ■ . 

All PAs in A' also appear in A, and all PA-indices of ff ' are PA-indices of ff . It follows 
that C',7' is also an elementary cycle of W,A, unless one of the pieces, say W! n ,,js 
contains the part of ff ' corresponding to H y q. More precisely, this happens when n(i{) = 
y and «'(//) = «(0). (In Figure 7, «'(28) = 1 = «(0) but «(28) = 32 =y.) But in that case, 
we may replace ff/ n , n ^ with W, ; V oH y xt o Wo.«(o)' an d choose the appropriate rotation 
(starting with 0), which yields an elementary cycle C,I of W,A with C w C'. This proves 
the first statement. 

The above construction maps elementary cycles of W',A' to elementary cycles of 
W , A. It is easy to see that they are all mapped to different elementary cycles, which only 
contain pieces that are subwalks of Wby, and that every elementary cycle of W,A that 
contains a piece of Wq y is covered this way. It remains to consider elementary cycles 
of W,A that contain a piece of W y Q. Since {Q,y} G A and n(y) = 0, there is only one 
such elementary cycle C,I (with C = W" , as defined above). This proves the second 
statement. □ 

As a first application of Lemma 19, we can determine the number of elementary 
cycles. 

Corollary 20 Let A be a proper assignment set for closed walk W in 3$, and let p be the 
number of PA-indices ofW. Then the number of elementary cycles ofW,A is p/2+h 

Proof: If A = then there are no PA-indices, and W itself is the only elementary cycle, 
which proves the statement. 

Otherwise we can use induction: choose a PA {x,y} G A with x = n(y). Such a 
PA exists since A is a non-crossing perfect matching on the PA-indices. W.l.o.g we 
may assume that x = 0, since considering a different rotation of W and changing A 
accordingly does not change p or the number of elementary cycles. 

Now consider W\ A' and ff" as defined above using {0,y}. Let p' denote the num- 
ber of PA-indices of ff'. ff" is an elementary cycle of ff, hence a right-turn walk, 
which has no PA-indices. Therefore every PA-index of ff other than or y corresponds 



to a PA-index in W', so p' + 2 = p. By Lemma 19 the number n of elementary cy- 
cles in W is equal to n' + 1, where n' is the numbers of elementary cycles of W',A'. 
Since A' is a proper assignment set for W' (Proposition 18) we may use induction, so 
n = n' + l=p'/2 + 2 = p/2+l. □ 

The above corollary will now be used to deduce that all elementary cycles defined 
by a proper assignment set have turning number 1, and hence admit a clockwise solution 
by Proposition 17. 

Lemma 21 Let A be a proper assignment set for a closed walk W in 33. Then t(W) = 1 
if and only if every elementary cycle C ofA,W has t(C) = 1 

Proof: The number of PA-indices of W is denoted by p. The number of left turns 
that W makes at non-PA-indices is denoted by Left* (W). Since W makes a left turn at 
every PA-index, we have 

6t(W) = Right(W) - Left(W) = Right(W) -p - Left*(W). 

Let c £ denote the set of all elementary cycles of A,W. When summing the difference 
between right and left turns over all elementary cycles we obtain 

£ (Right(C) - Left(C)) = Right(W) - Left*(W) + 2p. 

Here we used the following observations, (i) Elementary cycles make right turns at PA- 
indices of W, and every PA-index contributes a right turn to two elementary cycles, (ii) 
Non-PA-indices of W contribute the same type of turn to one elementary cycle, (iii) 
All indices of elementary walks that do not correspond to pieces of W correspond to 
internal vertices of assignment paths; these vertices of assignment paths contribute a 
left turn to one elementary cycle and a right turn to another, hence these terms cancel. 
Combining this with | = p/2 + 1 (Corollary 20) yields 

6t(W) + 3p = Right(W) - Left*(W) + 2p = 

£ (Right(C) - Left(C)) = £ (Right (C) - Left(C) - 6) + 3p + 6 <s=> 

t(W)= £( f (c)-l) + l. 

Because t(C) — 1 > for all Cef (Proposition 16, Proposition 17), this proves the 
statement. □ 



The proof of the next lemma is illustrated in Figure 8 (using the same example as 
before). 

Lemma 22 Let W, y, W' and W" be as defined above. If S' and S" are clockwise 
solutions of W' and W" respectively, then a clockwise solution S of W exists with 
£?(S)=jrf(S') + {0,y} 



Fig. 8. Combining two partial solutions. 



Proof: Let S' = (H' ,B') be a clockwise solution to W = W , y oH yfi , and S" = 
(H",<j)",B") be a clockwise solution to W" = W y floH 0j . Let x = k — y (recall that k 
is the length of W), so W^' x = W yfi . 

Let Z?', = Vq,Vj, . . . ,v(, and let B" Q = Vq, v", . . . , v". Note that these paths indeed 
have the same length / and 

f(vO = 0"(vL)forj = O,...,/. (1) 

Since {0,y} is a PA of W, we have I > 1. and y are PA-indices, so W makes left 
turns at and y, and therefore W' makes right turns at and y, and W" makes right 
turns at and x. Because both S' and S" are clockwise solutions, d(v' ) = d (v{) = 2 and 
c/(vq) = d(v',') = 2 follows. Since these two paths are mapped to the same paths in 28 
but in reverse direction, and both 5' and S" are clockwise solutions, we see that v\ has 
degree 2 when v 1 !. has degree 3 and vice versa, for i = 1, — 1 . 

Construct H by starting with a copy of H' and a copy of H", and for all i £ {0, . . . , /}, 
identify the vertex of H' with the vertex v" '_■ of H". Call the resulting vertex v,. For 
i G {0, ...,/ — 1}, this results in two parallel edges between v; and v,-+i. Delete one 
edge of every such parallel pair. Because I > 1, this gives again a 2-connected graph. 
Using the above observations on the degrees of v\ and v", we see that d(vi) = 3 for all i: 
either a vertex of degree 2 is identified with a vertex of degree 3 and two incident edges 
are removed (when i = 1 ,...,/— 1), or two vertices of degree 2 are identified and one 
incident edge is removed (when i = 0, /). 

Choose an embedding of H in which every inner face of H' or H" is again an inner 
face of H, and which has boundary cycle vi,B' 1: ... ,B' y l ,\'o,B", . . . ,B"_ 1: vi. A LIH 
from H is constructed by setting (j)(u) = (j>'(u) for all u & V(H) C\V{H') and (j>(u) = 
<j)"(u) for all u e V(H)nV(H"), and 0(v,-) = 0(v-) = ^(vf_!) for the new vertices (see 
Equation (1)). 

We observe that H is a hexagonal patch: all inner faces of H correspond to inner 
faces of H' or H", and thus all have length 6. The degree constraints hold since the new 
vertices v; have d(vi) = 3 and boundary vertices of H' and H" that were not identified 
remain boundary vertices. 



Next we show that <j> is aLIH. For vertices of H that have no neighbor in V (H')\{v' , . 
or no neighbor in V(H")\{v'q, ... , v'/}, the local injectivity follows from the local injec- 
tivity of (j)' and <j)" . The only two vertices of H for which this does not hold are vo and 
v/. But the three neighbors of V; in H are mapped to W\, Wj_i and a vertex on Ho y , 
which are all different. A similar statement holds for vq. This proves that is again a 
LIH. 

We observed above that the boundary cycle of H is 

B = vi,B 1 ,.. . ,B y _ l ,vo,B 1 , . . . ,E x _ x ,Vi, 

which is mapped by exactly to Wo,Wi, ... ,W y -i,Wy,W y+ i, ... ,Wfc_i,Wo = W. So 
(H,(j),B) is a solution to W. Since we started with clockwise solutions, it is also obvi- 
ously a clockwise solution. (By Lemma 1 we do not have to check this for all boundary 
vertices.) 

An assignment path from Z?J to B'j in H' yields an assignment path from B, to Bj in 
H. In addition, the path vo, . . . , V/ in H is mapped by to H v ,o, and no internal vertices 
of this path lie on the boundary of H. Hence this solution assigns to y. Now all PA- 
indices of W are accounted for (W" contains no PA-indices), which shows that s/ (S) = 
^(S') + {0,y}. □ 

Finally, we have collected all the ingredients that are necessary to prove the remain- 
ing lemma. 

Proof of Lemma 9: We show that if a proper assignment set A is given for a walk 
W with t(W) = 1, then a clockwise solution S of W with s^(S) =A exists. 

If A = then since A is a proper assignment set, this implies that W has no PA- 
indices and W itself is the only elementary cycle of A, W. This is a right-turn walk by 
Proposition 16, and since t(W) = 1, it has a clockwise solution S by Proposition 17. 
Note that £/(S) = = A since there are no PA-indices. 

If A 7^ 0, then choose an arbitrary {x,y} G A with n(y) = x (this exists since A is a 
non-crossing perfect matching). We may again assume w.l.o.g. that x = 0. Use this to 
define W', W" and A' as before. A' is a proper assignment set for W' (Proposition 18). 
Lemma 19 shows that for every elementary cycle C' of W',A', an elementary cycle C of 
W,A with C«C exists. By Lemma 21, t(C) = 1 follows from t(W) = 1, so t(C) = 1 
for every such elementary cycle C'. Therefore by applying Lemma 21 again, we obtain 
t(W') = 1. At this point we have a new closed walk W' in £% with turning number 1, 
with fewer PA-indices, and a proper assignment set A' for it. Hence by induction, W' 
admits a clockwise solution S' with £/(S') = A'. W" itself is an elementary cycle of 
W, so t(W") = 1 (Lemma 21), and therefore W" admits a clockwise solution S" as well 
(Proposition 16, Proposition 17). Now Lemma 22 shows that S" and S" can be combined 
into a clockwise solution S for W with £/ (S) =A. □ 



C Proofs of Section 5 

Proof of Lemma 13: We show that subroutine S(i,j) correctly calculates Su when the 
stated S Xi y values are known, in time 0(d(i)). Clearly, the algorithm only uses values 
S x , y with — 1 <y~ x < j — i for the calculations. 



Observe that throughout the algorithm, the value of m equals the size of some PM 
of Gy. We now show that in some line, the size of a MPM is considered, which proves 
correctness. 

Let M be a MPM of Gy. If j < i then M = which is considered in line 1. Now 
suppose j > i. If M contains no edge incident with i, it is a PM of G,+i ; and considered 
in line 2. Otherwise, let iv € M. If v = j, then M consists of this edge and a PM of 
Gj+i t—\, which is considered in line 6. Otherwise, because M is non-crossing, it can be 
partitioned into a MPM of Gy and G v+ i ,-, which is considered in the for-loop. 

The complexity of the algorithm is determined by the for-loop, which iterates at 
most d(i) times. □ 

Proof of Lemma 14: We show that subroutine N(i,j) correctly calculates iVy when 
the stated N Xl y and S x ,y values are known, in time 0(d(i)). The algorithm only uses 
values S Xi y with — 1 <y — x < j — i and values N x . y with — 1 < y — x < j — i for the 
calculations. 

Clearly, line 1 returns the correct answer (in this case the empty set is the unique 
MPM). Otherwise, the algorithm adds x to the number N whenever x different MPMs 
of Gy are found that have not been considered earlier. We show that all cases are 
considered and no MPMs are double counted, which shows that the correct answer is 
returned in line 7. 

If Sij = Sj + \j then all MPMs of G,-+i ; are also MPMs of Gy, which explains 
line 3. This accounts for all MPMs of Gy that do not contain an edge incident with i. 
Gy admits MPMs M that contain an edge ij if and only if Sy = +5/+i j-i. Then 
M — ij is a MPM of G,+i. 7 _i, and there is a bijection between such MPMs, which 
explains line 6. It remains to consider MPMs M of Gy which contain an edge iv with 
2 + 1 < v < j — 1. These can be decomposed into a PM of G l v and a PM of G v +ij, 
which must be MPMs, so Sy = Sy + S v +i. 7 -. In fact, if Sy = Sy +S v +ij then every 
combination of a MPM of Gy and a MPM of G v+ i j gives a unique MPM of Gy that 
contains iv, which explains line 5. All cases are now considered, so at the end of the 
algorithm, N = Nij. 

The complexity of the algorithm is determined by the for-loop, which iterates at 
most d(i) times. □ 



